Coexistence of gapless excitations and commensurate charge-density wave in the 

2H-transition metal dichalcogenides 
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An unexpected feature common to 2H-transition metal dichalcogenides (2H-TMDs) is revealed 
with first-principles Wannier functions analysis of the electronic structure of the prototype 2H-TaSe2: 
The low-energy Ta "5d,2" bands governing the physics of charge-density wave (CDW) is dominated 
by hopping between next-nearest neighbors. With this motivation we develop a minimal effective 
model for the CDW formation, in which the unusual form of the hopping leads to an approximate 
decoupling of the three sublattices. In the CDW phase one sublattice remains undistorted, leaving 
the bands associated with it ungapped everywhere in the Fermi surface, resolving the long-standing 
puzzle of coexistence of gapless excitations and commensurate CDW in the 2H-TMDs. 
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Charge-density waves (CDWs) in solids has been a 
topic of central interest in condensed matter physics for 
many years Q. Recent scanning tunneling microscopy 
experiments showing a periodic modulation in the local 
density of states in cuprate superconductors 0, Q , has 
reinvigorated such interest. Despite being one of the ear- 
liest discovered class of materials which exhibit a CDW 
at low temperatures many properties of the 2H- 
transition metal dichalcogenides (2H-TMDs) are still not 
understood, leading to much recent theoretical 0, H, 
and experimental jl, H [H 111^1 III 13 13 research 
effort (for a review, see Ref. [lal)- Two key issues con- 
cerning the CDW phase in these materials deserve the 
most attention. First, controversy exists between dif- 
ferent experimental groups on the driving mechanism of 
CDW originating from quantitative differences between 
the angle-resolved photoemission spectroscopy (ARPES) 
data. While some experimental results 
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suggest that the hexagonal Fermi surfaces around the 
r point are consistent with the CDW nesting vector, 
others 0, Q, 0| indicate that this Fermi surface is 
too large to give the correct nesting vector. Second, 
and of a more qualitative nature, ARPES measurements 
[1 [Hill [Hill El [13 find no evidence of a gap opening 
on the hexagonal Fermi surface, in direct contrast with 
traditional wisdom of CDW materials. 

In this Letter, we focus on the latter issue and sug- 
gest a simple picture for why such gapless excitations are 
permitted in the CDW phase. Using density-functional 
theory, the electronic structure of prototype 2H-TMD, 
2H-TaSe2, is analyzed with a newly developed Wannier 
function (WF) approach [l7l[li|. and a striking feature 
is revealed: the low-energy bands near Fermi surfaces, 
which governs the physics of CDW, is dominated by 
hopping between second- nearest neighbors. This special 
nature of hopping, in combination with the triangular 
lattice vectors, effectively splits the system into three 
weakly coupled triangular sublattices. Since the CDW 



state gives distortion of only two of the sublattices, whose 
stability is illustrated with a simple model, such unique 
electronic structure naturally leaves the bands associated 
with the undistorted sublattice ungapped in this CDW 
phase, resolving the puzzle of the observed gapless exci- 
tations along the nested regions of the Fermi surface. 

The lattice structure of 2II-TMDs consists of stacked 
layers of 2D-triangular lattices of transition metals (e.g. 
Ta or Nb) sandwiched between layers of chalcogen atoms 
(e.g. S or Se). Rough estimation of the ionization leads to 
the Ta^+Se^^ configuration with one valence bd^i elec- 
tron left per Ta atom that forms the metallic bands at the 
Fermi level. Indeed, the band structure from our first- 
principles calculations |E3 (see Fig. ^ shows a strong 
character in the two metallic bands corresponding to 
two weakly coupled TaSe2 sandwiches per unit cell. Un- 
expectedly, little Se p characters are found in these two 
bands. The calculated low-energy bands agree well with 
experiments, except that the saddle bands on the FK and 
AH lines are not as flat and close to the Fermi level as 




FIG. 1: First-principles band structure (dots) with the 
(black circles) and dxy/d^2_y2 (blue circles) characters shown. 
The bands below —0.7 eV are mainly Se p bands. Also shown 
are the bands constructed from low-energy WFs (green solid 
lines) and a 2D 'nesting' model (red dotted line; see text). 
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FIG. 2: a) Low-energy WF centered at Ta sites with ag (d^^) 
symmetry, colored to show its gradient from positive (red) to 
negative (blue). Notice the local dxy / d^2 _y2 symmetry in the 
hybridization tail (circled) located at neighboring Ta sites, 
b) Schematics of the WF in the layer of the Ta triangular 
lattice, with a similar color scheme giving the sign of the WF 
and three small dots marking the positions of Se atoms in 
the next layer, c) and d) Schematics of phase interference in 
hopping to first and second nearest neighbors, respectively. 
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FIG. 3; (a) The Fermi surface from the tight-binding band 
structure. White indicates unoccupied states and black in- 
dicates occupied states, (b) The Fermi surface for a slightly 
smaller chemical potential showing the extended saddle bands 
(small white regions). The first Brillouin zone is the hexago- 
nal cell formed by the red dashed lines. 
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Based on the first-principles ground state, the low- 
energy Hilbert space can be accurately extracted via lo- 
cal WFs (see Fig.E)), which we constructed by extending 
recently developed energy- resolved method llal to ii^- 
corporate desired symmetry As expected, the WF 
located at each Ta site has strong (1^2 symmetry near the 
center, before extending its unusual tails of dxy / d^^ _y2 
symmetry to the nearest neighboring Ta sites due to 
strong hybridization with the dxy/dx2_y2 orbitals near 
the K and H points (see Fig. This particular shape 
of the WF results in an intriguing feature in the hop- 
ping integral (evaluated via = ^/i u"^ with den- 
sity functional theory Hamiltonian, h^^^ , and Wannier 
states l/i) and \v)). That is, the second neighbor hop- 
ping, <2 = 115 mcV, overwhelms the first neighbor hop- 
ping, ti = 38 meV, due to remarkable phase cancelation 
in the latter case (to be discussed in more detail below). 
In addition, interlayer hoppings are found to be compa- 
rable to first-neighbor in-plane hopping with t^.i = 29 
mcV and tj_^2 = 23 meV. 

A simple microscopic picture for the unexpected dom- 
inance of second-neighbor hopping can be obtained from 
the symmetry of the WFs. As shown in Fig. \^c) and 
(d) , the contributions to the hopping parameters between 
neighboring WFs come mainly from overlap of their hy- 
bridization tails, since the tail-center [dz2 — dxy/dx'^-y'^) 
overlap gives negligible contribution due to its odd parity. 
While the first-neighbor hopping suffers seriously from 
the phase cancelation (illustrated by the opposite sign in 
Fig. |2Ic)), the second-neighbor hopping benefits greatly 
from the phase coherence (the same sign in Fig. |2Id)) of 
the overlap. Such symmetry consideration should hold 
for all 2H-TMDs of the same class, due to their similar 
local environment around the transition metal sites. 

Specifically in 2H-TaSe2, this unusual electronic struc- 



ture provides a plausible intuitive resolution to the puz- 
zling experimental observation of gapless excitations in 
the CDW phase. Indeed, with the dominating second 
neighbor hopping in a triangular lattice, the system ef- 
fectively splits into three weakly coupled sublattices. As 
we discuss below, one of the sublattices remains undis- 
torted in the CDW phase (see Fig. 0)) and therefore the 
bands associated with it are ungapped. 

With the dominance of the second neighbor hopping 
established, we now move on to construct a minimal, 
'nesting' model which captures the essential physics of 
the gapless CDW in 2II-TMDs. We start with the sim- 
ple 2D tight-binding energies given by 

£°k = EV|COs(k.R) (1) 

R 

where R runs over the triangular lattice defined by lat- 
tice vectors ai = a(\/3/2, l/2)and a2 = a(V3/2, -1/2). 
In addition to t2, te = ^2/3 [23| (all other hoppings ne- 
glected) is introduced to produce the Fermi surface of an 
almost perfectly nested "hexagonal checkerboard" pat- 
tern similar to the recent ARPES data with extended 
saddle bands (very close to but below the Fermi energy 
for extended regions along FK) |i[ll[Il[l3, as shown in 
Fig. El The corresponding band (with t^ adjusted to 140 
mcV) compares reasonably well with the first principles 
results (see Fig. ^ the two metallic bands are degenerate 
in the 2D model). As we show below, even with such a 
perfect Fermi surface nesting, no gap is opened upon the 
formation of CDW. 

Continuing the development of our minimal model, we 
next consider the CDW lattice distortions. The detailed 
neutron diffraction experiments of Moncton et al. |2^ 
have determined that the ionic displacements have Ei 
symmetry, which corresponds to longitudinal motion of 
the Ta atoms in the basal plane with amplitude given 
by experiment. However, the fitting procedure to the 
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measured geometric structure factors was insensitive to 
the overall phase of the distortions. Following this 
work, the atomic displacements having Si symmetry cor- 
responding to the triple period CDW in the 2H-TMDs 
are given by 



(5R 



u cos(Q • R - 

Q 



(2) 



Here u is the amplitude and the sum runs over the vec- 
tors Qi = bi/3, Q2 = b2/3, and Q3 = -(bi + b2)/3, 
where the reciprocal lattice vectors are given by bi = 
^(1/V3, 1) and ba = ^(1/^3,-1). The above atomic 
displacements also splits the lattice into three indepen- 
dent sublattices, where one of these sublattices does not 
experience displacements for any (p. We will determine 
the overall phase factor tp by minimizing the total energy. 
It can be seen that the magnitude of the displacements 
given by Eq. 10) will not depend on Lp. Thus the elas- 
tic energy of the system will not depend on the phase 
of the CDW for this model system and our problem is 
reduced to finding the phase that minimizes the energy 
of the conduction band. Expanding the crystal poten- 
tial to first order in (5R given by Eq. Q leads to the 
perturbation 



k,Q 
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h.c. 



(3) 



where the wave vector k is summed over the first Bril- 
louin zone. For simplicity, we assume that the change in 
the hopping parameters due to the lattice distortion is 
proportional to the change in the absolute distance be- 
tween neighboring atoms: (J^RR' oc ((5R — (5R') • (R — R'). 
Then 
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l)e-^''-«-Q • R (4) 



where the lattice vector R is summed over the sec- 
ond nearest neighbors to the atom at the origin and 
7|R|=2 > is the clcctron-phonon coupling constant in 
the unit of energy/distance. Since the unit cell of the 
distorted lattice contains nine sites of the original lattice, 
the renormalized energies ej^"^ are given by the eigenval- 
ues of a 9 X 9 matrix which is rather cumbersome and we 
will give its explicit expression elsewhere. 

We can now write the total energy as a function of the 
amplitude and phase of the distortion as 



(5) 



where p(e) = X^k ?i '^(^ ~ ^L"^) density of elec- 

tronic states and Ec\{u) is the elastic energy which, as 
we noted before, is independent of the phase of the dis- 
tortion. The k-intcgration is performed by using a fine 
Monkhorst-Pack mesh and corresponding weights [3| in 
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FIG. 4: Atomic displacement pattern corresponding to the 
phase </5 = 7r/2 in Eq. (|2J. The inset shows the total energy 
of the CDW state as the function of phase tp where the total 
energy units are arbitrary and the phase is in units of tt. 



the irreducible Brillouin zone of the triangular lattice. 
We determine the chemical potential /i by considering 
two extreme cases: (i) Fixed particle number iV, which 
corresponds to an isolated metallic band at the Fermi en- 
ergy, (ii) Fixed = 0, which corresponds to significant 
spectral weight from the other bands at the Fermi en- 
ergy. The realistic situation, with two metallic bands at 
the Fermi energy arising from the multi-layer structure 
should reside somewhere between these two extremes. 
We find that the unanimous minimum for both extreme 
cases occurs when p = 7r/2. This minimum will be- 
come more pronounced with increased electron-phonon 
coupling constant. In Fig. ^ we show the corresponding 
atomic displacement pattern, which is consistent with the 
charge maxima seen in scanning tunneling microscopy ex- 
periments [13, ■ Furthermore, to check the robustness 
of this result, we have performed the same calculation, 
but with tg set to zero, and have found that the minimal 
total energy still occurs at </? = 7r/2. 

Now that the phase of the CDW has been determined 
we will analyze the renormalized quasiparticle disper- 
sion in the presence of the CDW. The energy spectrum 
from the undistorted sublattice is shown in Fig. EJa) 
and that from distorted sublattices is shown in Fig. El^b) 
along the FM direction, which is along the nested re- 
gion of the F-centered hexagonal Fermi surface. Clearly, 
those associated with the undistorted sublattice do not 
change [Fig. El^a)]. Thus, the corresponding sub-bands 
remain metallic in the CDW phase. On the other hand, 
the bands originating from the two distorted sublat- 
tices are doubly degenerate and we find a gap opens 
at the Fermi energy [Fig. IHb)]. Moreover, to make 
comparison with experiment more direct, in Fig. Etc) 
and Fig. EJd) we present theoretical ARPES spectra 
A(k, w) = i/(ct;)ImG(k, Lj) where G'(k, w) is the single 
particle Green's function and f{uj) is the Fermi distribu- 
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FIG. 5: The bands in the low temperature CDW state orig- 
inating from the undistorted (a) and distorted (b) sublat- 
tices. The increase in the number of bands corresponds to 
backfolding resulting from the 3x3 supercell. Theoretical 
ARPES spectra for the normal (c) (corresponding to vanish- 
ing electron-phonon couphng) and CDW (d) states for wave 
vectors over the nested region of the Fermi energy along FM. 
The energy range is —0.47 to eV relative to the Fermi energy 



and the momentum range is 



to 



tion function for wave vectors in a small region of the 
nested portion of the Fermi surface. To emulate exper- 
imental data, we have chosen a broadening of 77 = 40 
mcV of the spectral density function. In Fig. E^b) both 
the gapped and ungapped bands arc visible (for an ac- 
count of the weightings of satelite bands in CDW mate- 
rials, see Voit et al. |27l|'). The most direct comparison 
between experimental data is with the work of Valla et 
al. ^3 where there is a plot similar to Fig. El^a). It is 
also shown in this paper that no gap opens along FK, 
being consistent with Fig.|3Ib). 

A question that naturally follows is, how robust this 
result is when the finite first neighbor interaction — which 
mixes all bands and thus destroys the exact decoupling 
into the three independent sublattices — is taken into ac- 
count. We examine this issue, switching on the first- 
neighbor electron-phonon coupling constant 71 in Eq. Q), 
up to a third of the value of 72 as suggested by the first 
principles results: ti ~ ^2/3. As expected, we find that 
the degeneracy of the bands originating from the dis- 
torted sublattices shown in Fig. E^b) is lifted. In addi- 
tion, the triple degeneracy originating from the undis- 
torted sublatticc shown in Fig. Efa) at the Fermi energy 
is lifted. However, this does not produce a quasiparticle 
gap at the Fermi energy. More specifically, this triple 
degeneracy is lifted in such a way that two of the en- 
ergies are increased to above the Fermi energy and the 
other one is decreased to below the Fermi energy. It can 
be seen that this will indeed not gap the Fermi surface, 
given 71 is considerably smaller than 72. 



In conclusion, we have studied the CDW state in the 
2H-transition metal dichalcogenides and found that due 
to a unique feature in the electronic structure of these 
materials revealed from first-principles calculations, the 
triangular lattice can be effectively decoupled into three 
independent sublattices, with one remaining undistorted 
in the CDW phase. As illustrated with a model calcu- 
lation, this leads to the remarkable situation where no 
regions of the entire Fermi surface become gapped even 
when these materials exhibit a commensurate CDW. 
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